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1. Introduction 

Nuclear shadowing (NS) in deep-inelastic scattering (DIS) off nuclei is usually studied via 
nuclear structure functions. Assuming shadowing region of small Bjorken xsj < 0.01 the structure 
function F2 per nucleon turns out to be smaller in nuclei than in a free nucleon. This fact affects 
then the corresponding study of nuclear effects mainly in connection with the interpretation of the 
results from hadron-nucleus and heavy ion experiments. 

There are different treatments of NS depending on the reference frame. In the infinite mo- 
mentum frame of the nucleus NS can be interpreted as a result of partem fusion [jlj fz|] leading to a 
reduction of the parton density at low Bjorken xbj- In the rest frame of the nucleus, however, this 
phenomenon looks like NS of the virtual photon hadronic fluctuations and is occurred due to their 
multiple scattering inside the target (see refs. p, |], and references therein, for example). 

In the rest frame of the nucleus, the dynamics of NS is controlled by the time scale known 
as the effect of quantum coherence. It results from destructive interference of the amplitudes for 
which the interaction takes place on different bound nucleons. It can be estimated by relying on 
the uncertainty principle and Lorentz time dilation. For the lowest Fock component of the photon 
it reads, 

2v 

where v and Q 2 is the photon energy and virtuality, respectively and M qq is the effective mass of 
the qq pair. It is usually called coherence time, but we also will use the term coherence length (CL), 
since light-cone (LC) kinematics is assumed, l c = t c . CL is related to the longitudinal momentum 
transfer q c = \/l c . 

Note, that higher Fock states containing gluons, \q~qG), \qq2G), ... have larger effective masses 
than qq fluctuation resulting in shorter coherence times. 

The Green function approach ^ is a very effective tool for the study of NS naturally incor- 
porating the effects of CL. Such a formalism is based on the solution of the evolution equation for 
the Green function. Usually, one has been forced to obtain an analytical solution leading to a har- 
monic oscillator form (see Eq. ( |2. 1 1| )). It requires, however, to implement several approximations 



into a rigorous quantum-mechanical approach like a constant nuclear density function (2.19) and a 



specific quadratic form (2.18) of the dipole cross section. 

We remove such approximations for the \qq) Fock state solving the evolution equation numer- 
ically using algorithm from ref. [JJ]. It is supported also by the fact that the predictions for NS [||] 
corresponding to the lowest Fock component of the photon showed quite large differences from 
approximate calculations [Q, ^ using the harmonic oscillator Green function approach for l c < 
(Ra is the nuclear radius). 

Calculations of NS presented so far in the color dipole approach [f7|, ||, were performed 
assuming only qq fluctuations of the photon and neglecting higher Fock components containing 
gluons and sea quarks. The effects of higher Fock states are included in energy dependence of 
the dipole cross section, Oq q (r,s) 1 . However, investigating nuclear effects, these Fock states lead 
also to gluon shadowing (GS), which has been neglected so far when the model predictions were 
compared with experimental data. 

'Here r represents the transverse separation of the qq photon fluctuation and s is the center of mass energy squared. 
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It was demonstrated in ref. [Q] that a contribution of GS to overall NS becomes effective at 
small xbj < 0.01. It confirms the necessity for inclusion of GS in the model predictions in the 
shadowing region ~ 0.0001 < x$j < 0.01, where the data exist. 



2. Light-cone dipole approach to nuclear shadowing 



In the rest frame of the nucleus NS in the total virtual photoabsorption cross section af^ (xBj,Q 2 ) 
(or in the structure function F^ixBjjQ 2 )) can be decomposed over different Fock components of 
the virtual photon. Then the total photoabsorption cross section on a nucleus can be formally 
represented in the form 



o£t A (x BJ , Q 1 ) = A <t/ / v (x B j, Q 2 ) - Ao tot (x Bj , Q 2 



jfN 



(2.1) 



where Ao tot (x B j, Q 2 ) = Aa tot (qq) +Ao tot (qqG) + Ao tot {qqlG) + x B j is the Bjorken variable and 
°iof ' ( x Bj, Q 2 ) is total photoabsorption cross section on a nucleon 



oL N (,xbj, Q 2 ) = j d 2 r J da ¥^(r, a, Q 2 ) a qq {f, s 



(2.2) 



The dipole cross section Oq q (r,s) in Eq. (2.2) represents the interaction of a qq dipole of trans- 
verse separation r with a nucleon [§J. It depends on the cm. energy squared s. The dependence 
of Gq q (r,s) on r and energy is universal and flavor independent allowing to describe various high 
energy processes in an uniform way. <jq q (r,s) is known to vanish quadratically oc r 2 as r —> due 
to color screening (property of color transparency [||, |9|]). It cannot be predicted reliably because 
of poorly known higher order perturbative QCD (pQCD) corrections and nonperturbative effects. 
However, Oq q if,s) can be extracted from experimental data on DIS and structure functions using 
reasonable parametrizations. In this case pQCD corrections and nonperturbative effects are natu- 
rally included. 



There are two popular parametrizations of Oq q (r,s), GBW presented in [10] and KST sug- 
gested in Detailed discussion and comparison of these two parametrizations can be found in 
refs. (Tl], ^]. Whereas GBW parametrization can not be applied in the nonperturbative region of 
small Q 2 , KST parametrization presents a good description of the transition down to limit of real 
photoproduction, Q 2 = 0. Because we study the shadowing region of small xbj < 0.01 and avail- 
able experimental data from the E665 and NMC collaborations cover small and moderate values of 
Q 2 < 2 4- 3 GeV 2 we prefer the latter parametrization. 

The KST parametrization [FJ] has the following form containing an explicit energy dependence, 



Gqqfas) = a oO) 



exp 



r 2 



*8M 



where 



o»(s) = o? ?(s)[\ + 



3g(£) 



(2.3) 



(2.4) 



In Eq. (2.3) the energy-dependent radius Rq(s) 



1000GeV z . In Eq. (O) <J*f(s) = 23.6 (s/s 



\0.079 



0.88 (s/soy^^fm with X = 0.28 and ^ = 
- 1.432 (s / sq)~ 0A5 mb corresponding to the 
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Pomeron and Reggeon parts of the %p total cross section [12] and (r 2 h ) K 



0.44 fm 2 represents the 



mean pion charge radius squared. 

Several comments about KST parametrization are in order: 



i) The KST dipole cross section (23) is proportional to r 2 for r — > 0, but flattens off at large r. 

ii) The saturation scale Rq{s) decreases with increasing s. The energy dependence of o qq {r,s) 
correlates with r. At small r the dipole cross section rises with a hard pomeron intercept, 0.36, and 
at large separations it still depends with a soft pomeron intercept, 0.08, on energy. 



iii) In contrast to GBW parametrization [10], energy dependent ansatz, Eq. (2.3), allows to 
describe simultaneously hadron-hadron scattering and DIS. This improvement at large separations 
leads to a worse description of the short-distance part of the dipole cross section which is responsi- 
ble for the behavior of F^XBj^Q 2 ) at large Q 2 . To satisfy Bjorken scaling the dipole cross section 
at small r must be a function of the product sr which is not the case for the parametrization in 



Eq. (£3|). 

iv) The form of Eq. ( [23] ) successfully describes the data for DIS at small xbj only up to 
Q 2 lOGeV 2 . Nevertheless, this interval of Q 2 is sufficient for the purpose of the present paper. 



The second ingredient, m qq (r, a, Q ), in (22) is the perturbative distribution amplitude ("wave 
function") of the qq Fock component of the photon. It depends on the photon virtuality Q 2 and 
the relative share a of the photon momentum carried by the quark. The corresponding form for 
transversely (T) and longitudinally (L) polarized photons can be found in refs. [|T|, |14L Ejh: 



<<f(r,a,<2 2 ) 



271 



(2.5) 



where % an d X are the spinors of the quark and antiquark, respectively; Z q is the quark charge, 
Nc = 3 is the number of colors. Ko(er) is a modified Bessel function with 



a(l 



a)Q 2 + m 2 



(2.6) 



where m q is the quark mass. The form of the operators T L can be found in 



The distribution amplitude Eq. (2.5) controls the transverse qq separation with the mean value, 



1 



2 2 a(l-a) + m 2 



(2.7) 



For very asymmetric qq pairs with a or (1 — a) <m 2 /Q the mean transverse separation 
(r) ~ l/m q becomes huge since one must use current quark masses within pQCD. A popular recipe 
to fix this problem is to introduce an effective quark mass m e ff ~ Aqcd which should represent the 
nonperturbative interaction effects between q and q. However, we introduce this interaction explic- 
itly using corresponding phenomenology based on the LC Green function approach developed in 



The Green function G qq {r2,z.2',r\,z\) describes the propagation of an interacting qq pair be- 
tween points with longitudinal coordinates z\ and z.i and with initial and final separations f\ and r\. 
This Green function satisfies the two-dimensional Schrodinger equation, 



• d r 
l lz~^ 



,(r 2 ,z 2 ;n,zi) 



2va(l -a) 



+ V qq {z2,r 2 ,OC) 



G qq (r 2 ,Z2;n,zi) , 



(2.8) 
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with the boundary condition 



Gc iq (r 2 ,z2;n,zi)\ Z2 = zl = d (n-r 2 



(2.9) 



In Eq. (2.8) V is the photon energy and the Laplacian A r acts on the coordinate r. 

Studying propagation of a qq in vacuum the LC potential V qq (z2,r2,cc) in (2.8) contains only 
the real part, which is responsible for the interaction between the q and q. For the oscillator form 
of this potential, 



ReVsJz2,f 2 ,a) 



a 4 (a)r 2 2 



(2.10) 



2va(l-a) 

one can obtain the analytical solution of the evolution equation ( j2~8| ) represented by the harmonic 
oscillator Green function [15], 



Gq q (r 2 ,Z2\r\,Z\) 



where Az = z 2 — Z\ and 



a 2 (a) 



2 71 i sin(ft) Az 



exp 



ia 2 (a) 
sin(ft)Az) . 



v \ +r 2 )cos(co Az.) — 2f\-r 2 

ie 2 Az 

xexp 

[ 2va(l-a) 



CO 



a 2 (a) 



(2.11) 



(2.12) 



v a(l - a) 

The shape of the function a(a) in Eq. ( 2.10| ) is presented and discussed in ref. [Hj. 

Matrix element (2.2) contains the LC wave function squared, which has the following form 
for T and L polarizations in the limit of vanishing interaction between q and q, 



%(f,a,Q 2 ) 



2N c a e , 

m 2 



^ q (r,(X,Q 2 ) 



£ Zj [m 2 f K (e,r) 2 + [a 2 + (1 - a) 2 } z 2 K x (e r) 



/=i 



(2nf 



£ Z 2 f Q 2 a 2 {\-a) 2 K Q (ei 



(2.13) 



(2.14) 



where and K\ are the modified Bessel functions. 

The effects of higher Fock states \qq), \qqG), \qq2G), ... are implicitly incorporated into 
the energy dependence of the dipole cross section Oq q (r,s) naturally included in realistic KST 
parametrization Eq. (|Z3|). 

Investigating DIS on nuclear targets, the corresponding formula for NS, Aa tot (xBj,Q 2 ) = 
Ao to t{qq), representing the shadowing correction for the lowest qq Fock state, has the following 
form [0] 

Ae tot (x Bj ,Q 2 ) = -Re J d 2 b J dzi Pa(P,zi) J dz 2 p A (k>,z 2 ) da A(z u z 2 ,a) , (2.15) 



with 



A(zi,Z2,a) = J d 2 r 2 »F*- 9 (r 2 ,a,e 2 ) O qq (r 2 ,s) J d 2 r x G qq (r 2 ,z 2 ;n,zi) O qq {r u s) ^ qq (r u a , Q 2 ) . 

(2.16) 
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In Eq. ( |2.15| ) p A {b,z) represents the nuclear density function defined at the point with longitudinal 



coordinate z and impact parameter b. 

The shadowing term A.o tot (x B j,Q 2 ) = Aa tot (qq) in ( |2.15| ) is illustrated in Fig.[I[ At the point z\ 



the initial photon diffractively produces the qq pair (y*N — > qqN) with transverse separation f\ . The 
qq pair then propagates through the nucleus along arbitrary curved trajectories, which are summed 
over, and arrives at the point zi with transverse separation r^. The initial and final separations 
are controlled by the LC wave function of the qq Fock component of the photon ^ 9 (r, a, Q 2 ). 
During propagation through the nucleus the qq pair interacts with bound nucleons via the dipole 
cross section Oq q {r,s), which depends on the local transverse separation r. The Green function 
Gqq{r2,Z2\r\ ,z\) describes the propagation of the qq pair from z\ to zi- 

Describing propagation of the qq pair in nuclear medium the Green function Gq q {r2,Z2\r\,z\) 



satisfies again the time-dependent two-dimensional Schrodinger equation (2.8). However, the po- 
tential in this case acquires in addition also the imaginary part responsible for attenuation of the qq 
photon fluctuation in the medium and has the following form 

ImVqq(z 2 ,7,a) = - ° qq [ r,s) p A (b,z 2 ) • (2.17) 



harmonic oscillator potential Vq q (r) oc r 2 , i.e 



As we already mentioned above the analytical solution of Eq. ( |2.8| ) is known only for the 



Oqq{r )S )=C(s)r 2 , (2.18) 

and uniform nuclear density 



p A (b,z)=p &(R 2 A -b z -z z ) (2.19) 



in Eq. (2.17) for the imaginary part of the LC potential. 
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In this case the solution of Eq. (fO|) has the same form as Eq. (2.1 1), expect that one should 
replace oo CI and a 2 (a) ==> b(a), where 



b(a) ^/a A (a)-ip A {b,z)va(\ — a)C{s) 



(2.20) 



va(l — a) va(l-a) 
Determination of the energy dependent factor C(s) in Eq. ( [2. IS ) and the mean nuclear density 



Po in Eq. ( |2.19| ) can be realized by the procedure described in \q, g]. 

In the general case when there is no restrictions for l c , if l c ~ Ra one has to take into account 
the variation of the transverse size r during propagation of the qq pair through the nucleus. The 
overall total photoabsorption cross section on a nucleus o^ A = Oj A + e' ol[ A , assuming that the 
photon polarization e' = 1 . If one takes into account only qq Fock component of the photon the 
full expression after summation over all flavors, colors, helicities and spin states has the following 



form [16] 



° fA (xBj,Q Z 



A^ N (x Bj ,Q 2 )-Aa(x Bj ,Q 2 ) 







2 






1 d 2 r dao qq (r,s) 1 




+ 


%j(r,a,Q 2 ) 


') 



N c a, 

W) 2 f- 



Nf 1 

~V r r°° poo p 1 p p 

J^Zjfte J d 2 b J dzi J dz 2 J da J d 2 r x J d 2 r 2 



:p A {b,z\) PA(b,z 2 ) a qq (r 2 ,s)a qq (r u s)< a 2 + (l-a) 



r\r 2 



+ 



m 2 +4Q 2 a 2 (l-a) 2 K (e n) K Q (e r 2 ) } G qq {r 2 ,z 2 \?uZi) . (2.21) 



*F^l L (r, a, Q 1 ) is given by Eqs. (2.13) and (2.14), respectively. 



■ m 



The shape of 

In the high energy limit, l c ^> Ra, the transverse separation r between q and q does not vary 
during propagation through the nucleus (Lorentz time dilation). Consequently, the kinetic term in 



Eq. (2.8) can be neglected and the Green function reads 

G qq (b\r 2l z 2 \r u z\)\v^°° = 8(r 2 -ri) exp 



1 f Z2 
~O qq (r 2 ,s) / dz,p A (b,z) 



(2.22) 



After substitution of the expression (2.22) into Eq. (2.21), one arrives at the following results: 
a,JWe 2 ) = 2 / d 2 b I d 2 r I da\l-ex V ^~(j qq (r,s)T A (b) 



o 



2N c a e , 
{2Kf 



-£z|| a 2 + (l-a) 2 K x {erf + m 2 f + AQ 2 a 2 (1 - a) 2 K 2 (er) j, (2.23) 



where 



T A (b)= dzp A (b,z) 



(2.24) 



7 



Nuclear shadowing in the light-cone dipole approach 



Jan Nemchik 



is the nuclear thickness calculated with the realistic Wood-Saxon form of the nuclear density with 
parameters taken from [fn|]. 

At photon polarization parameter e' = 1 the structure function ratio F^Xbj, Q 2 ) /^(xsj, Q 2 ) 
can be expressed via ratio of the total photoabsorption cross sections 

Fj(xBj,Q 2 ) 4 A (xBj,Q 2 ) + o[ A (x B j,Q 2 ) 

F 2 N (x S j,Q 2 ) aJf N (x B j,Q^) + ar N (x B j,Q2) J ^ } 

where the numerator on right-hand side (r.h.s.) is given by Eq. ( [2.2 1[ ), whereas denominator can be 



expressed as the first term of Eq. (2.21) divided by the mass number A. 

The nonperturbative q — q interaction is included by replacements Ko(er)/2n =>- <J>o(£,r,A) 
and Ki(er)/2n =>■ <S>i(e,r,A) in all perturbative expressions. The corresponding functions <&o,l 
and parameter A are defined in ref. 



We solve the evolution equation for the Green function, Eq. ( |2.8| ), numerically using algorithm 
from ref. Such an exact solution is performed for realistic KST parametrization of the dipole 
cross section ( |2~3| ) and nuclear density function in the realistic Wood-Saxon form with parameters 
taken from ref. [|l7|]. 

Finally we would like to emphasize (see also the next section) that the qq Fock component 
of the photon represents the higher twist shadowing correction [Q], and vanishes at large quark 
masses as l/mj. This does not happen for higher Fock states containing gluons, which lead to 
GS. Therefore GS represents the leading twist shadowing correction [§J. Moreover, a steep energy 
dependence of the dipole cross section Cq q {r,s) (see Eq. (^3|)) especially at smaller dipole sizes r 
causes a steep energy rise of both corrections. 



3. Gluon shadowing 

Investigating nuclear effects at small xbj < 0.01, the higher Fock states of the photon contain- 
ing gluons also contribute to NS. Because of a shorter coherence time (lifetime) of these fluctua- 
tions, GS will be dominated at higher energies, i.e. at smaller xbj- values than NS coming from 
the lowest \qq > Fock component. However, no data for GS are available and one should rely on 
calculations. 

NS for \qq) Fock component of the photon is dominated by T photon polarizations, because the 
corresponding photoabsorption cross section is scanned at larger dipole sizes than for the L photon 



polarization. The transverse qq separation is controlled by the distribution amplitude Eq. (2.5 ), with 



the mean value given by Eq. ( J2.7| ). Contributions of large size dipoles come from the asymmetric 
qq fluctuations of the virtual photon, when the quark and antiquark in the photon carry a very large 
(a — ► 1) and a very small fraction (a — > 0) of the photon momentum, and vice versa. The LC wave 



function for L photons (2.14) contains a term a (1 — a) , which makes considerably smaller the 



contribution from asymmetric qq configurations than for T photons (see Eq. ( |2.13 )). Consequently, 



in contrast to T photons, all qq dipoles from L photons have a size squared x l/Q and the double- 
scattering term vanishes as oc l/Q 4 . The leading-twist contribution for the shadowing of L photons 
arises from the \qqG) Fock component of the photon because the gluon can propagate relatively 
far from the qq pair, although the q-q separation is of the order l/Q 2 . After radiation of the gluon 



8 



Nuclear shadowing in the light-cone dipole approach 



Jan Nemchik 



the pair is in an octet state and consequently the \qqG) state represents a GG dipole. Then the 
corresponding correction to the L cross section is just GS. 

Interpretation of GS is reference frame dependent. In the infinite momentum frame this phe- 
nomenon looks analogical as gluon-gluon fusion. Within a parton model interpretation, the gluon 
clouds of nucleons which have the same impact parameter overlap at small x B j in longitudinal di- 
rection. This allows gluons originated from different nucleons to fuse leading to a gluon density 
which is not proportional to the density of nucleons any more. This is GS. 

In the rest frame of the nucleus phenomenon of GS corresponds to the process of gluon radia- 
tion and shadowing corrections related to multiple interactions of the radiated gluons in the nuclear 
medium. This is a coherence phenomenon known as the Landau-Pomeranchuk effect, namely the 
suppression of bremsstrahlung by interference of radiation from different scattering centers. It de- 
mands a sufficiently long coherence time of radiation, a condition equivalent to demanding a small 
Bjorken x B j in the parton model. 

There are still very few numerical evaluations of GS in the literature, all done in the rest frame 
of the nucleus. GS can be identified as the shadowing correction to the L cross section coming 
from the GG dipole representing \q~qG) Fock component of the photon. For evaluation of GS is 
important to know about the transverse size of this GG dipole. This size has been extracted in 
ref. [Q] from data for diffractive excitation of the incident hadrons to the states of large mass, the 
so called triple-Pomeron region. It was found that the mean dipole size of the GG system (radius 



of propagation of the LC gluons) is rather small , ro ~ 0.3 fm J18Q. It results in a not very strong 
onset of GS. 

The smallness of the size of quark-gluon fluctuations has been incorporated via a nonperturba- 
tive LC potential in the Schrodinger equation for the Green function describing the propagation of 
a quark-gluon system. The strength of the potential was fixed by data on high mass (M|) diffrac- 
tion pp — > pX [§]. This approach allows to extend the methods of pQCD to the region of small 
Q 2 ^ Go = 4/ r o- At higher Q 2 > Qq GS slowly (logarithmically) decreases in accordance with 
expectations based on the evolution equation 

We repeated the calculations [Q] of the ratio of the gluon densities in nuclei and nucleon, 

„ , G A {x Bj ,Q 2 ) AOtot(qqG) 

Rg[XB - Q ] ~ AG N (x Bj ,Q 2 ) ~ 1 ' (3J) 

where Ao to t (qqG) is the inelastic correction to the total cross section of ^ related to the creation of 
a \qqG) intermediate Fock state. Further calculation details can be found in As an illustration 
of not very strong onset of GS, here we present Rg{x B j, Q 2 ), Eq. ([Q|), for different nuclear thick- 



nesses ?a(&). Using an approximation of constant nuclear density (see Eq. (2.19)), T^b) = PqL 



where L = 2^JR 2 i — b 2 , the ratio Rg(x B j, Q 2 ) is also implicitly a function of L. An example for the 
calculated L-dependence of Rc{ x Bj,Q 2 ) at Q 2 = 4GeV 2 is depicted in Fig. ^|for different values 

Of XBj. 

We calculated GS only for the lowest Fock component containing just one LC gluon. Inclusion 
of higher multigluon Fock components is still a challenge. However, their effect can be essentially 
taken into account by eikonalization of the calculated Rg{ x bj,Q 2 ), i-C- the dipole cross section, 
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Figure 2: The ratio of the nucleus-to-nucleon gluon densities as function of the thickness of the nu- 
cleus, L — T(b) /po, at Q 2 = 4 GeV 2 and different fixed values of xsj- Figure is taken from ref. Jl~ij ] . 



which is proportional to the gluon density at small separations, should be renormalized everywhere, 



RgOs 



(3.2) 



According to Eq. ( p.2[ ) we will demonstrate that GS suppresses the total photoabsorption cross 
section on a nucleus o^t (xBj,Q 2 )- We expect a non-negligible effect of GS in the shadowing 
region of small xbj < (0.001 -j- 0.01 ) and at small and medium values of Q 2 ~ 2 -f - 3 GeV 2 corre- 
sponding to the kinematic range of available data. 



4. Numerical results 

Here we present the available data vs. realistic predictions for NS in DIS based on exact 
numerical solutions of the evolution equation for the Green function corresponding to the lowest 
qq Fock component of the photon. Such a comparison is performed for the shadowing region of 
small XBj < 0.01. 

We take into account also a contribution of GS which should increase overall nuclear sup- 
pression. The effects of GS are calculated for the lowest Fock component containing just one LC 
gluon. Inclusion of higher Fock components with more gluons is realized by eikonalization of the 



calculated Rg( x bj,Q ), i-e using renormalization ( pT2| ) 



For a numerical solution of the Schrodinger equation for the Green function we adopt an 
algorithm from ref. []5]]. Because available data from the E665 Jl9| ] and NMC |2C] collaborations 
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-l 



Figure 3: Nuclear shadowing for lead. Calculations correspond to exact numerical solution of 
the evolution equation for the Green function, Gq CI (r2,Z2',r\,zi), for the lowest Fock component of 
the photon using KST Q parametrization of the dipole cross section and realistic nuclear density 
function of the Woods-Saxon form [p"7|]. The solid and dashed curves represent the predictions 
calculated with and without contribution of gluon shadowing, respectively. 



cover the region of small and medium values of Q 2 < 4 GeV we use the KST parametrization of 
the dipole cross section [§] which is valid down to the limit of real photoproduction. 

Performing a numerical solution of the evolution equation for the Green function, the imagi- 
nary part of the LC potential ( 2. 17[ ) contains the corresponding KST dipole cross section as well. 
The nuclear density function Pa(^,z) is taken in the realistic Wood-Saxon form with parameters 
from ref. [|l7[|. Because available data from the E665 and NMC collaborations correspond to very 
small values of Q 2 < 1 GeV 2 at small xsj < 0.004, the nonperturbative interaction effects between 
q and q are included explicitly via the real part of the LC potential of the form ( 2.10 ). 

Effects of NS are studied via xbj- behavior of the ratio ( 2.25[ ) divided by the mass number A. 
Firstly we present NS for the lead target in Fig. || at different fixed values of Q 2 . The solid and 
dashed curves represent the predictions obtained with and without contribution of GS, respectively. 
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One can see that as a manifestation of a shorter l c for higher Fock states, the onset of GS 
happens at smaller xgj than the quark shadowing. Fig. || clearly demonstrates not very strong onset 
of GS in the range of x^j £ (0.01,0.0001) where the most of available data exist. Besides the 
effects of GS are stronger at smaller Q 2 because corresponding Fock fluctuations of the photon 
have a larger transverse size. 
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Figure 4: Comparison of the model with experimental data from the E665 [ |l9| | and NMC [ |20| l 
collaborations. Calculations correspond to exact numerical solution of the evolution equation for 
the Green function, Gq q (fi-,Zi\f\,z\), for the lowest Fock component of the photon using KST [§ 
parametrization of the dipole cross section and realistic nuclear density function of the Woods- 
Saxon form [jl7||. The solid and dashed curves are calculated with and without contribution of gluon 
shadowing, respectively. 



Saturation of NS at low xgj < 10 at the level given by Eq. (2.23) is realized only for energy 
independent dipole cross section (see parametrization (2.18)). However, it is not so for the realistic 
energy-dependent KST parametrization, Eq. (2.3). 

In Fig. [| we present a comparison of the model predictions with data from the E665 [ |l^ ] and 
NMC [20] collaborations. Fig. shows a quite reasonable agreement with experimental data in spite 
of absence of any free parameters in the model. One can see that the effect of GS produces an 
additional NS, which rises with the mass number A. Fig. ^| also demonstrates a non-negligibility of 
GS already in the region of xbj < 0.001 -j- 0.01 . Very large error bars especially at small xsj ~ 10 
do not allow to investigate separately the effect of GS. Therefore more accurate new data on NS in 



12 



Nuclear shadowing in the light-cone dipole approach 



Jan Nemchik 



DIS at small xgj are very important for the further exploratory study of nuclear modification of the 
structure functions and also GS. 
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Figure 5: Comparison of the model predictions for the ratio Pb/nucleon obtained without (upper 
thick solid line) and with GS (lower thick solid line) with other models, versus xgj at fixed Q 2 = 
3 GeV 2 . HKM are the results from FS-NLO from [g, Bartels from @, Sarcevic from 
Frankfurt from [Ell] (Q 2 = 4 GeV 2 ), Armesto from & EKS98 from 0] and ACKLS from | 



Finally, Fig. || represents a comparison of the model predictions for NS with the results from 
other models at Q 2 = 3 GeV 2 (except the results of ref. [21] which are at Q 2 = 4 GeV 2 ). At small 
Xbj = 10~ 5 we predict quite large effects of GS (compare upper and lower thick solid lines). Note 
that the difference between models rises towards small values of XBj following from a different 
treatment of various nuclear effects. 

The models presented in Fig. |5] can be divided into several groups: 

i) models where parton distribution functions (PDF's) has been determined from data using 
Dokshitzer-Gribov-Lipatov-Altareli-Parisi (DGLAP) evolution equation 

Here the nuclear PDF's have been parametrized performing next to leading order (NLO) [22] or 
leading order (LO) [ 23 , 24 ] global analysis of nuclear DIS and Drell-Yan data. 

ii) models based on Glauber-like rescatterings 

Here the model in ref. [25 ] is based on an application of a saturating ansatz for the total /"-nucleon 
cross section in the proton. This ansatz is then extended to the nuclear case by its introduction 
in a Glauber expression. Within the model from [^] a Glauber ansatz provides with the initial 
condition for DGLAP evolution. 

iii) models based on Gribov inelastic shadowing 

Here in ref. [E7l] nuclear structure functions are studied using relation with diffraction on nucleons 
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which arises from Gribov's reggeon calculus. The model presented in [21] employs again some 
parametrization of hard diffraction at the scale <2q, which gives nuclear shadowing through Gribov's 
reggeon calculus similarly as in ref. [ 27]. Then nuclear suppression calculated at Qq is used as 
initial condition for DGLAP evolution. 

iv) models based on high-density QCD 
Here the model in ref. [ p8| ] is based on numerical solution of a non-linear equation for small-^gj 
evolution and on application of this equation for the case of nuclear targets. 

Fig. |5] shows large differences in predictions of NS at small xbj among different models. It 
gives a stimulation to obtain new more accurate data on nuclear structure functions by lepton-ion 
collider planned at RHIC [E^]. It can help us to discriminate among different models. 



5. Summary and conclusions 

We present a short review of the color dipole approach based on the LC QCD Green func- 
tion formalism, which naturally incorporates the interference effects of CT and CL. Within this 
approach [0, |], ||j] we study NS in DIS at small Bjorken xbj- 

Calculations of NS corresponding to qq component of the virtual photon are based on an 
exact numerical solution of the evolution equation for the Green function. It allows to use realistic 



parametrizations of the dipole cross section (GBW [ 10] and KST M) and realistic nuclear density 



function of the Woods-Saxon form [ 17] 



Because available data from the shadowing region of xbj < 0.01 coming mostly from the 
E665 and NMC collaborations cover only small and medium values of Q 2 < 4 GeV 2 , we prefer 
KST parametrization [Q] of the dipole cross section. On the other hand the data obtained at a lower 
part of the xbj- kinematic interval correspond to very low values of Q 2 < 1 GeV 2 (nonperturbative 
region). For this reason we include explicitly the nonperturbative interaction effects between q 



and q taking into account the real part of the LC potential Vq q ( |2.10| ) in the time-dependent two 



dimensional Schrodinger equation ( |2.8| ) 



In order to compare the realistic calculations with data on NS, the effects of GS are taken into 
account. The same path integral technique [||] is applied in this case. GS was calculated only for 
the lowest Fock component, \qqG). Effect of higher Fock components containing more gluons was 
essentially taken into account by eikonalization of the calculated Rg(xbj, Q 2 ) using renormalization 



(3.2). We demonstrate that the onset of GS starts to be effective at xbj ~ 0.01. It rises towards 
small xbj because higher Fock components with more gluons having shorter coherence time will 
contribute to overall NS. Such a situation is illustrated in Fig. |3[ 

Performing numerical calculations, we find that our model is in reasonable agreement with 
existing experimental data (see Fig. Q). Large error bars and incompatibility of the experimental 
results from the E665 and NMC collaborations do not allow to study separately the effect of GS. 
Therefore more accurate new data on NS in DIS off nuclei at still smaller xbj < 10~ 5 are very 
important for the further exploratory study of GS effects. 

Comparison among various models shows large differences for the Pb/nucleon ratio of the 
structure functions at XBj = 10~ 5 and Q 2 = 3 GeV 2 (see Fig. Q). It has a large impact on the cal- 
culation of high-pr particle spectra in nuclear collisions at RHIC and LHC. Such large differences 
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at small xgj among different models should be testable by the new more precise data on nuclear 
structure functions which can be obtained by lepton-ion collider planned at RHIC 
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